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The Standard Model has been successful in describing many fundamental aspects of particle physics. However, there are some 
remaining puzzles that are not explained within the context of its present framework. We discuss the possibility to discover 
new physics in the ATLAS Detector via a four-fermion contact interaction, much in the same way Fermi first described Weak 
interactions. Using a simple ratio method on dimuon events, we can set a 95% C.L. lower limit on the effective scale A =7.5 TeV 
(8.7 TeV) for the constructive Left-left Isoscalar Model of quark compositeness with Jz? = 100 pb~' (200 pb~') of data at 
y/s = 10 TeV. 



1. Introduction 

New physics discovery may be just around the corner as 
the Large Hadron Collider (LHC) prepares for first colli- 
sions this winter. The Standard Model has thus far shown 
impressive predictive power and agreement with experi- 
ment; yet the cause of Electroweak Symmetry Breaking 
(EWSB), necessary to give mass to the and Z° bosons, 
remains experimentally unconfirmed. In the last decades, 
many models outside of the current Standard Model frame- 
work have been developed to address this. 

One way to model a new interaction with unknown cou- 
plings is in the form of a "4-fermion" contact interaction. 
Interactions with a dimuon final state have been chosen for 
this analysis, as they provide a clean signature in the early 
stages of understanding the ATLAS detector. 

In the first year of LHC running, we expect 100 - 
200 pb~' of data at y/s = 10 TeV. In these proceedings, 
we discuss the feasibility of discovering or setting limits 
on new physics via contact interactions with early ATLAS 
data. Note also that the results presented here have not yet 
been officially approved by the ATLAS Collaboration. 



2. Theoretical Background 

Beyond the Standard Model processes, such as large ex- 
tra spacial dimmensions (LED) in the ADD model QJJ] or 
quark/lepton compositeness [2], may be described as a 4- 
fermion contact interaction. In the same spirit as the Fermi 
Interaction describes j3 -decay without directly knowing 
the intermediate process [3], one can write an effective La- 
grangian describing a new interaction [0]: 

+ lLR{WL7 fl VfL){Y'R7nWR) 

+ riRLiwf^i^tiW'Ll ] , (1) 

where g is a coupling constant, and r and y/' L R are the 
incoming and outgoing left and right fermionic fields, re- 



spectively. The interaction appears experimentally as a de- 
viation from the SM dilepton mass spectrum, which origi- 
nates from Drell- Yan (DY) production (gg— > 7,Z^ l + l~). 
The value of rj is the sign of the interaction term; there can 
be either constructive (tj = —1) or destructive (rj = +1) 
interference with the DY process. 




(a) (b) 

Figure 1: Representative Feynman diagrams for various 4- 
fermion contact interactions. Examples of possible production 
at the LHC are shown for (a) dijet production and (b) dilepton 
production. 

Some examples of contact interaction diagrams are 
shown in Fig. Q] In the case of Fig. QJ, the "contact" 
can be understood as being between incoming partons and 
final-state charged leptons at some scale A. More specifi- 
cally, the Lagrangian for the 4-fermion contact interaction 
with a dimuon final state (qqjj.jj.) is given by: 

g 2 

% = -^\nu,{qLfqh)\JiL7^L) 

+ r\RR{m7 ll m){V-L7ixV-L) 
+ f]LR (qLT^tfL) {p-r7h^r) 
+ riRL(qR7 >1 qR){P-L7^L} } • (2) 

Here, qi r are the left- or right-handed quark doublets, u^r 
and cIl,r are the left or right quark singlets, and /1l,r are the 
left- or right-handed muon singlets. 

The total cross section is then given by the general form 
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Figure 2: Production mechanism of Drell-Yan with additional 
contact term with scale A in the dimuon final state. 



Table I Cross-sections for benchmark A values in the construc- 
tive Left-left Isoscalar Model for pp collisions at \fs = 10 TeV 
(for Mft^ > 120 GeV). 



A (TeV) 


a x bf [pb] 


5 


13.28 


7 


12.85 


9 


12.75 


12 


12.54 


oo(DY) 


12.52 



with a DY term (the cross-section of the Standard Model 
Drell-Yan spectrum), an interference term / which goes 
as A~ 2 , and the contribution of the contact interaction C 
which goes as A~ 4 . 



3. Contact Interaction Analysis Method 

Monte Carlo event generation using PYTHIA J5| has 
been done for the parity-conserving Left-left Isoscalar 
Model (L-LIM) of fermion compositeness, corresponding 
to the first term in Eq. ©. We chose four benchmark sig- 
nal values to calculate the expected limit: A = 5, 7, 9 and 
12 TeV. In the event generation stage, the following selec- 
tion criteria were used: 

• two or more final state muons from the hard scatter- 
ing 

• dimuon invariant mass M^ > 120 GeV 

• transverse momentum pj > 5 GeV and pseudo- 
rapidity 1 7] | < 2.8 for final state muons. 

These requirements were chosen in order to increase statis- 
tics in the signal region and to be within the geometri- 
cal acceptance of the Muon Spectrometer. In Table U we 
show the production cross-section times the muon branch- 
ing fraction (X — > jiji) for each of the benchmark A val- 
ues. Detector response was simulated using the standard 
ATLAS fast simulation, and the resulting dimuon invariant 
mass was calculated. We then multiply by a mass depen- 
dent k-factor (ranging from 1.31 to 1.15 from 120 GeV to 
2000 GeV [6]) and a dimuon reconstruction efficiency fac- 
tor of (0.85) 2 = 0.07225. After event generation and fast 
simulation, muon candidates are chosen to be within pseu- 
dorapidity |tj| < 2.5 (the geometrical acceptance of the 
Inner Detector [7]) and dimuons are required to have an 
invariant mass > 120 GeV. Fig. [3] shows the result- 

ing dimuon differential cross-section in the constructive 
Left-left Isoscalar Model. The DY spectrum corresponds 

to A — * oo. 
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Figure 3: Dimuon invariant mass spectra for various A values. 
Note that the signal distribution tends to the DY shape as A — > °°, 



above an invariant mass threshold (Mo) and divide by the 
number of events in a control region where the expected 
shape and cross-section of the DY are tightly constrained 
from previous experiments: 



R 



KDY,d — 



Nm>m 
Nm<m 



(4) 



Note that Ra-^™ = Rdy, the expected ratio from the DY in- 
variant mass spectrum. The number of events in any region 
is given by N = Jz?eobf(X — * jUjll), where J2? is the in- 
tegrated luminosity, e is the acceptance and reconstruction 
efficiency, a is the differential cross-section from Eq (O, 
and (bf(X — > fifJ.)) is the branching fraction in the dimuon 
channel. When taking the ratio Ra.dy from Monte Carlo 
simulations or Rj from collected data, the uncertainty in 
the luminosity cancels in the ratio. The effect of strongly- 
correlated uncertainties between the two regions is also re- 
duced, which is especially advantageous in the first months 
of data taking, when the uncertainties are expected to be 
large. 

The ratio in Eq. (@}, expected for the DY-only spectrum, 
can be compared with an expected signal ratio using the 
following significance calculation: 



3.1. Limit Setting and Discovery 
Potential Using the Ratio Method 

One method to look for an excess of events in the high- 
mass DY spectrum is to simply count the number of events 



S/i, 



Ra — Rdy 



°1 + °dy 



(5) 



Here ox is the statistical uncertainty on the signal ratio, 
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assuming that new physics did exist but had fluctuated 
downwards and was missed. Systematic uncertainties on 
the background ratio are incorporated in Ody, which in- 
clude detector effects and our limited knowledge of Stan- 
dard Model parton density functions, backgrounds, etc. In 
calculating the actual limit, we would use real data in place 
of the DY-only Monte Carlo sample for comparison to the 
ratio expected for a given A scale. The A scale which re- 
sults in a significance of Sn,„ = 1.96 corresponds to a 95% 
confidence level limit. 

3.2. Expected Reach With Early Data at 
v^= 10TeV 

The effect of systematic uncertainties on the signifi- 
cance calculation was determined for the DY background. 
It was found that muon resolution, muon efficiency, 
and k-factor uncertainties were the largest contributors 
to the overall ratio uncertainty, as they are momentum- 
dependent. The efficiency uncertainty was chosen to 
rise from 1% to 15% over the mass range of 120 GeV- 
2000 Ge V, the resolution uncertainty calculation was based 
on a mis-alignment the Muon Spectrometer four times 
worse than design, and the k-factor uncertainty was based 
on minimizing and maximizing the slope of the mass- 
dependent k-factor functional form within errors previ- 
ously determined |g|. Note that the values quoted are an 
exaggeration of expected first-year systematics. Table UTI 
summarizes these results. 

Table II Expected systematics for the first months of ATLAS 
running. N^igh and N{ ow are the expected number of events for 
100 pb~' above and below Mo, respectively. 



Mass cut @ 450 GeV 


Ni ow 


Nhi g h 


DY ratio 


% diff 


Nominal SM 


660.5 


7.6 


0.0115 




Resolution 


732.7 


9.0 


0.0123 


6.7% 


Pt scale (1% up) 


699.3 


7.9 


0.0113 


-1.5% 


pr scale (1% down) 


632.5 


7.3 


0.0116 


1.1% 


Efficiency (up to 15%) 


669.5 


7.8 


0.0117 


2.0% 


Efficiency (down 15%) 


651.4 


7.3 


0.0112 


-2.1% 


k-factor (max slope) 


635.3 


7.4 


0.0117 


1.9% 


k-factor (min slope) 


700.5 


7.8 


0.0111 


-3.1% 



To determine the impact of statistical uncertainty in 
the signal ratios (<7a), the RMS of Rj^ distributions from 
performing 10000 pseudo-experiments were calculated on 
each of the benchmark values (Fig. @1 Table [TTTil. The 
number of sample events (A/j) chosen for the i th pseudo- 
experiment was determined from a Poisson distribution 
about the expected number of events N expectec i in 100 pb~' 
(or 200 pb -1 ) of data. Each pseudo-experiment was then 
created by randomly sampling Nj events from a dataset cor- 
responding to hundreds of fb~ . 

The value of the expected ratios Rdy.a are highly depen- 
dent on the choice of Mo, the cutoff between the high and 
low dimuon invariant mass regimes. To account for this, 
we calculate Su m as a function of Mo (Fig. [5]). The Mo 
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Figure 4: Comparing various distributions for Jzf = 100 pb 
(top) and 200 pb~' (bottom). 



which maximizes Su m for each A value is used when set- 
ting the limit on A (see Table [TlTb - We have found that the 
maximum is very broad in Mo, and that increasing the lu- 
minosity does not affect the chosen mass cut significantly. 



Table III Ratios Rqy and Rj^ are shown for each A with corre- 
sponding best Mq cut in 100 pb . 



A 


Best Mq 


Drell-Yan only 


Drell-Yan + signal 


tC 


Rdy 






OX 


5 TeV 
7 TeV 
9 TeV 
12 TeV 


450 GeV 
540 GeV 
750 GeV 
780 GeV 


10.20 

805.52 
4.98 


0.0127 
0.0062 
0.0016 
0.0013 


44.33 

820.95 
13.66 


0.0541 
0.0170 
0.0046 
0.0025 


0.0082 
0.0046 
0.0024 
0.0017 


809.37 
1.27 


808.37 
3.73 


812.36 
1.07 


820.56 
2.00 


812.53 


805.93 



Note that in early data, while the resolution uncertainty 
is the largest systematic, statistical uncertainty dominates 
in calculating the ratio (Fig. [6). After fitting the signifi- 
cance for the various A values (5, 7, 9 and 12 TeV), we 
expect to be able to set a new limit on the constructive L- 
LIS model of fermion compositeness at A = 7.5 TeV with 
100 pb" 1 , or A =8.7 TeV with 200 pb" 1 (Fig.|7). 
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Figure 5: Distributions of significance (5/, m ) by varying Mo for different A values, with the blue curve corresponding to 100 pb and 
the red curve corresponding to 200 pb~ 1 . 




Lambda Figure 7: The expected 95% C.L. limit on A corresponds to a sig- 

nificance (Su,„) of 1.96o~, shown here for 100 pb -1 and 200 pb~' . 

Figure 6: Ratio versus Lambda, using Mq for maximized signif- 
icance. 
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